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Chapter 1 

Rappels 


Ce cours fait suite au cours de Mecanique des Fluides de l’annee L3-GSI. Ce premier chapitre vise donc ä rappeier quelques 
principes importants qui pourront etre necessaires par la suite. 

1.1 Proprietes des fluides et quantites physiques necessaires ä leur etude 

Definition d’un fluide : on appelle fluide un corps qui s’ecoule : liquide, gaz, rnais aussi des corps plus complexes tels que les 
polymeres ou les fluides alimentaires. Un fluide se caracterise par un certain nombre de grandeurs qui le caracterisent tels 
que sa densite, sa viscosite, etc 

• masse volumique : s’exprime en kg/m 3 . 

• densite : eile rnesure le rapport de la masse volumique du fluide rapportee ä celle de l’eau : sans unite. 

• viscosite cinematique et dynamique : la viscosite cinematique v est homogene ä une diffusivite ( m 2 /s ). La viscosite 
dynamique g est le produit de la precedente par la masse volumique:// = pv (en Pa.s ). 

• pression : se rnesure en Pascal. 

• Vitesse : en m/s 

• Acceleration : en m/s 2 

• Contrainte : homogene ä une pression (Pa ou N/m 2 )). 

1.2 Rappels d’hydrostatique 

On appelle hydrostatique un probleme caracterise par un fluide au repos. La vitesse est donc nulle en tout point soit :V = 0. 
Nous allons donc considerer un fluide au repos sournis ä un champ de force / par unite de volurne et exprimer la relation qui 
lie le champ de pression au champ de force. Apres Integration, cette relation locale donnera la valeur du champ de pression 
en tout point du fluide. 

On rappelle d’abord que pour un element de surface dS appartenant ä une surface fermee S entourant un volurne de controle, 
la force de pression que ressent dS est dirigee (par Convention) dans le sens oppose ä sa normale exterieure, soit: 

dfp = —p n.dS 

Ceci traduit le fait que les forces externes de pression ont tendance ä repousser la surface vers l’interieur du volurne qu’elle 
enferme. Pour une surface finie S , la resultante des forces de pression s’ecrira donc: 

— J J p ndS 

Considerons maintenant l’equilibre d’un volurne elementaire de fluide dx.dy.dz au repos et sournis ä un champ de force 
( fx'Jv'Jz ) (cf. Figure 1.1): 

La somme des forces auxquelles il est sournis vaut 0 de Sorte que l’on peut ecrire: 

/ — grad p = 0 
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Cette expression est homogene ä des N.m~ 3 et fournit l’occasion de rappeier la definition d’un gradient: 
le gradient est un Operateur qui s’applique ä une grandeur scalaire ( la temperature, la masse volumique, etc...) ; il traduit 
l’intensite des variations spatiales de cette grandeur. En coordonnees cartesiennes, il s’ecrit: 


grad, f 


dl 

dx 

d£ 

dy 

dl 

dz 


L’operateur grad se note aussi V : ainsi, Vp signifie grad p . 
En coordonnees cylindriques, l’operateur gradient est defini ainsi: 


grad f 


dl 

dr 

df 


dz 


Lorsque le champ de force / se limite ä la gravite (soit f z = — pg si on Oriente Taxe Oz vers le haut) on en deduit par 
projection selon Oz : —pg — ^ = 0 d’oü Ton tire: 


p + pgz = Cste 

Exemple 1: Force sur une paroi verticale 
Exemple 2 : Force sur une paroi inclinee 


1.3 Poussee d’Archimede 

C’est la poussee que ressent un corps lorsqu’il est immerge dans un fluide. Cette force intervient ä l’arret (bateau ou baigneur 
immobile) mais aussi en mouvement (une gouttelette d’eau qui tombe dans l’air est soumise ä son poids, aux forces de trainee 
mais aussi ä la poussee d’Archimede). 

Enonce : la poussee d’Archimede est egale au poids du volume fluide deplace. 

Lorsqu’on etliche une Situation dans laquelle la poussee d’Archimede intervient, il faut ecrire l’equilibre des forces (la 
somme des forces sur un corps vaut 0 que ce soit un corps immobile -cf. Exemple 1- ou un corps de vitesse constante c.a.d. 
d’acceleration nulle cf. Exemple 2). 

Exemple 1 : calcul de la position d’un iceberg d’eau douce dans l’eau salee 
Exemple 2: calcul de la vitesse limite de chute d’une spliere dans un fluide (viscosimetre ä chute) 


1.4 Rappels d’hydraulique 

Theoreme de Bernoulli : Ce theoreme exprime le fait que l’energie d’un fluide peut se transformer de 3 fagons : en pression, 
en energie potentielle ou encore en energie cinetique: 
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Si le fluide est parfait (cas ideal d’une viscosite nulle) ce theoreme exprime que la somme de ces 3 formes d’energie reste 
constante le long d’une ligne de courant , soit: 


P + ^pV 2 + pgz = Cste 

Une autre forme, obtenue apres division par pg pennet d’exprimer la Charge H : 

P V 2 

H = — -h — — “t - z — C ste 
P9 2 g 

qui s’exprime donc en metre. 

Dans la pratique, la viscosite du fluide est responsable de pertes d’energie de Sorte que la clrarge H s’epuise, et donc 
diminue, au für et ä mesure que le fluide avance dans les conduites. On distingue 2 types de pertes de Charge: 




Les pertes de cliarges regulieres qui traduisent la perte d’energie lorsque le fluide avance dans une conduite de diametre 
constant. Ces pertes sont donc proportionnelles ä la longueur de la conduite. On montre qu’clles s’ecrivent: 


D h 2g 


Dans cette expression: 

Dh est le diametre hydraulique defini par Dh = 7 ? ou Sest la section de passage du fluide et P le perimetre mouille. 
Vd est la vitesse debitante, rapport du debit volumique ä la section S. X est le coefficient de perte de Charge donne par 
les abaques (voir par exemplc l’abaque de Moody Figure 1.2). Ce coefficient n’est fonction que du nornbre de Reynolds 
Re = Vd t f h et de la rugosite relative k = oü e designe la hauteur (en m) des asperites. 

• Les pertes de Charge singulieres qui traduisent les pertes d’energie dans les accidents de canalisation (Te, coude, 
convergent, divergent, clapet, vanne, debitmetre, crepine, etc). Chaque perte de Charge est caracterisee par un coefficient 
K ä lire dans une table adaptee au type de perte de Charge considere. La perte de Charge s’exprime alors : 


V? 

A H sina = K— 
Slng 2g 


1.5 Rappels sur les machines hydrauliques 

Cette section vise ä rappeier la fagon dont on peut prevoir le debit d’une installation comportant une machine hydraulique 
(pompe ou turbine) donnee, ou , de fagon similaire, rappeier la fagon dont on choisit une machine hydraulique lorsqu’on veut 
faire circuler un certain debit fixe par le cahier des charges du pro j et par exemple. 

On ne considerera dans ce chapitre que deux types de machines hydrauliques: 


• Les pornpes : une pompe, comrne tout convertisseur d’energie, possede une caracteristique propre qui est donnee par 
le fournisseur. Cette caracteristique lie la Charge H fournie par la pompe au debit volumique Q qu’elle fera circuler. 
Pour la plupart des pornpes, et notannnent les pornpes centrifuges, la caracteristique H — Q est une courbe decroissante 
du fait que la Charge que peut vaincre la pompe diminue lorsque le debit augmente. Cette courbe est l’analogue de la 
courbe Tension-Courant d’une pile et a donc globalement la meine forme. 

Pour chaque pompe cette courbe depend des parametres suivants : 

- La vitesse de rotation de la pompe : lorsqu’on augmente la vitesse de rotation d’une pompe, le debit et la Charge 
augmentent de Sorte que la caracteristique H — Q se transforme en glissant vers la droite mais aussi vers le haut (Cf. 
figure 1.3) 

- Le diametre de la turbine : sur un meine corps de pompe , on peut monter plusieurs types de turbines de diametres 
differents. La figure 1.4 montre un exemple de caracteristique d’une pompe GEC ALSTOM KID P 250-525 : on y 
voit l’evolution de la caracteristique pour des diametres de roue allant de 434 mm ä 532 mm. Cette figure montre 
aussi l’evolution du rendement dont les isovaleurs forment ce que l’on appelle les collines de rendement Ce rendement 
est le apport de la puissance mecanique transmise au fluide sur la puissance mecanique regue par l’arbre de la pompe soit: 
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REGIME TURBULENT. TUYAUX RUGUEUX 


33 - Diagramme de Moody 
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NITOfeRO DE REYNOLDS R, 


Phyd QPljP- 

^ ~ ~P h _ /70 

1 arbre ^ ^ L 

Les courbes d’iso-puissance sur l’arbre figurent elles aussi sur la Figure mais sont en fait superflues puisse que la con- 
naissance de H 1 Q , et 77 suffit ä calculer la puissance. Enfin, la figure 1.4 indique la Charge minimale que l’on doit 
respecter ä l’aspiration de la pompe si on veut eviter la cavitation. Cette Charge se nomme le NPSH (pour Net Pressure 
Succion Head). Par exemple, avec de l’eau et pour un debit de 1000m 3 /h on lit un NPSH de 6m ce qui signifie que la 
pression ä l’aspiration ne doit pas etre inferieure ä 0,66ar. 

L’utilisation des pompes est tres repandue dans toutes les industries mais il est quand meine possible de classer leur 
utilisation en deux types: 

1: le fluide circule en circuit ferme (boucle de circulation pour chauffage, pour reaction chimique, filtration etc... 

2: le fluide est transvase d’un bassin ä un autre, generalement plus eleve : la pompe fonctionne en pompe de relevage 


• Les turbines : ce sont des recepteurs qui convertissent l’energie hydraulique en energie mecanique de rotation, alimen- 
tant ensuite un alternateur pour faire du courant ou, mais c’est maintenant rare, un moulin ou une scierie. Cette fois, 
plus le debit qui circule dans la turbine est grand et plus la Charge prelevee ä l’ecoulement sera grande : on a une 
courbe caracteristique H — Q croissante qui part de l’origine (cf. Figure 1.5). Dans une turbine, la vitesse de rotation 
n’est pas imposee comme dans une pompe : cela constitue une difference importante. La vitesse se calcule en ecrivant 
l’egalite entre la puissance fournie par l’arbre de la turbine et la puissance absorbee par la Charge de cette turbine , 
soit: 


Parbre — ( - - — TJturbine - PdQ P 


1 : 1 500 t/min 2: 2000 t/min 3: 2400 t/min 4: 2900 t/min 



Figure 1.3: Courbe caracteristique d’une pompe centrifuge :effet de la vitesse de rotation 


Calcul du point de fonctionnement d’une pompe Deux cas tres distincts se presentent: 

Cas 1: on connait le debit de la pompe car on l’a mesure (debitmetre ou autre methode): par simple lecture sur la car- 
acteristique de la pompe (voir par exemple figure 1.4) on en deduit la Charge de la pompe et son rendement, la puissance 
mecanique fournie au circuit hydraulique etc... 

Cas 2 : on ne connait pas le point de fonctionnement (cas le plus frequent) : on doit tracer la caracteristique du circuit et 
trouver graphiquement le point de fonctionnement. La caracteristique du circuit donne, ä partir des elements du circuits 
, et independamment de la pompe, la valeur de la perte de Charge (exprimee en metres) ä fournir en fonction du debit. 
Cette courbe est donc aussi une courbe H — Q dont l’intersection avec la caracteristique est le seul point de fonctionnement 
possible. 

Considerons la Figure 1.6 : la courbe 1 represente la caracteristique de la pompe. Imaginons qu’elle soit branchee sur un 
circuit ferme dont la caracteristique est donne par la courbe 2, globalement parabolique d’equation H = K.Q 2 : le point de 
fonctionnement sera le point A. Imaginons que l’on trouve le debit Qa trop grand : apres ajout d’une perte de Charge dans le 
circuit (fermeture plus marquee d’une vanne par exemple) la nouvelle caracteristique sera la coube 3 dont l’intersection avec 
la courbe 1 enfl donne un debit Qb inferieur ä Qa- Si le circuit de cette pompe est un circuit de relevage, la caracteristique 
sera teile que la courbe 4 qui pour Q = 0 represente la hauteur de fluide ä remonter et le point de fonctionnement sera en C . 
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Figure 1.4: Courbe caracteristique d’une pompe centrifuge ALSTOM 
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Figure 1.5: Courbe caracteristique d’une turbine 



Figure 1.6: Determination graphique du point de fonctionnement d’une pornpe 


1.6 Notions de cinematique 

Nous rappcllerons dans cette partie les notions de cinematique vues en L3. Choisissons une particule fluide comme entite 
elementaire (une particule fluide est un ensemble de molecules entourant un point M donne et supposees avoir la merne 
vitesse ä un instant donne). 

Nous allons voir que l’ecoulement peut etre decrit selon deux points de vue tres differents. 

Descriptions d’Euler et de Lagrange 

La description d’Euler consiste ä etablir ä un instant donne, l’ensemble des vitesses associees ä chaque point du fluide 
(Figure 1.7). 



Figure 1.7: Point de vue Eulerien : la vitesse, mesuree aux points M\ ou M% evolue au cours du temps. 


La vitesse Vm( 4) associee au point M evolue au cours du temps. A chaque instant l’ecoulement du fluide est decrit au 
moyen d’un charnp de vecteurs vitesse. Dans cette description, on appelle ligne de courant la courbe qui, en chacun de ses 
points, est tangente aux vecteurs vitesse (cf. Figure 1.8). 



Figure 1.8: Lignes de courant ä differents instants. 


La description de Lagrange consiste ä suivre une particule donnee au cours de son mouvement au sein du fluide. Ici, c’est 
l’evolution de la position des particules qui permet la description de l’ecoulement. On appelle trajectoire le lieu geometrique 
des positions successives occupees par une particule (cf. Figure 1.9). 

Lignes de courant et trajectoire sont bien deux notions tres differentes (cf. Figure 1.10), sauf si l’ecoulement est station- 
naire, c’est ä dire que le champ de vecteurs vitesse est constant dans le temps : il y a alors concidence entre lignes de courant 
et trajectoires. Nota : un ecoulement stationnaire est aussi denomme ecoulement permanent. 
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P(t 2 ) 


trajectoire de la 
,?f particule P 


P(t ± ) 


P(t 3 ) 


P(t 0 ) 


Figure 1.9: Definition de la trajectoire 


trajectoire 



Figure 1.10: Difference entre trajectoire et ligne de courant 
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Chapter 2 

Ecoulements de fluides parfaits 


L’etablissement des equations completes regissant l’ecoulement d’un fluide reel, c’est-ä -dire visqueux, est complexe. C’est 
pourquoi, dans ce chapitre, nous allons commencer par etablir les equations d’un fluide parfait ce qui suprime donc le terme 
lie aux contraintes visqueuses. Cette etape intermediaire permettra une transition entre la relation Fondamentale de la 
Dynamique (RFD) vue en Mecanique du Point et la Mecanique des Fluides oü on s’interesse ä ce qui se passe en un point 
donne de l’ecoulement. 


2.1 Definition d’un fluide parfait 

II s’agit donc d’un fluide tlieorique de viscosite nulle (/z = v = 0). On fait l’hypothese que les molecules du fluide glissent les 
unes sur les autres sans frottement (alors qu’il existe des forces d’interaction moleculaires dans un fluide reel). Signaions que 
l’Helium liquide ä quelques degres Kelvin se rapproche d’un fluide parfait tant sa viscosite tend vers 0. 

L’absence de viscosite se fait mieux ressentir lorsqu’on compare sur la Figure 2.1 les ecoulements de fluides reels et parfaits 
dans diverses situations. 


2.2 Rappel de la Relation Fondamentale de la Dynamique 

Pour une particule de masse m la RFD exprinre le lien qui existe entre la somnre des forces exterieures agissant sur m et son 
acceleration : ma = rn^r = f soit par unite de volume : 

pa = f v 

f v represente l’ensenrble des forces de volume qui peuvent agir sur un fluide : gravite, forces electriques, magnetiques, 
tlrermiques, etc 

2.3 Equation d’Euler 

La relation de la Dynamique permettrait d’etudier le nrouvenrent de clraque particule fluide une ä une mais le nornbre 
d’equations ä resoudre sera aussi grand que le nornbre de particules contenues dans le fluide. On realise donc une approclre 
dite de rnilieu continu dans laquelle on va exprimer la vitesse et l’accelration du fluide en un point donne de l’ecoulenrent 
(de coordonnees r et ceci en fonction du temps t). On clierche donc ä etablir une equation qui fournisse les valeur de la 
vitesse V(r,t) et de la pression P(f,t), ainsi que T(r, t), si le fluide n’est pas isotherme ou c(r, f) si le fluide transporte un 
constituant en solution de concentration C. 

Pour cela on considere donc un element de volume sur lequel on va d’abord exprimer la rsultante des foces de pression (cf. 
Figure 2.2. 

Pour une surface S et dans la seule direction x la force resultante (en Newton) vaut : 

dP 

+P{x)S — P(x + dx)S = — —Sdx 

dx 

soit en divisant par le volume fluide Sdx la force volumique selon Ox se resunre ä — En generalisant aux 3 directions de 
l’espace on deduit la force resultante de pression : 

fp = —gradP = — VP 


11 



Fluide parfait 


Fluide reel 


Ecoulement 
en conduite 



V=Cste pourtout r 
P=Cste pourtout x 




Laminaire 

V=V(r) 


Turbulent 


Ecoulement 
autour d'un cylindre 



a a 



a o o 


Trainee =0 


Trainee > 0 
tourbillons 


Ecoulement 
sur une paroi plane 


>/ 

► 


Trainee =0 


V=V(y) 
Trainee > 0 


Ecoulement 
dans un divergent 





AH=0 

pas de perte de Charge 


Recirculations 


Figure 2.1: Comparaison de l’ecoulement d’un fluide parfait et reel dans diverses situations. 


n 

P(x).S P(x+dx),S 

»n 





X + dx 


Figure 2.2: Bilan des forces de pression sur une tranche d’epaisseur x. 
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Nous allons maintenant exprimer l’acceleration. 
Expression Eulerienne de l’acceleration 





X + dx 


Figure 2.3: Ecoulement ID selon Ox. 


Considerons un ecoulement 1 D selon Ox : en M la vitesse vaut V x et en M' vaut V x . L’acceleration a x (suivant Ox) au 
sens de Lagrange, c’est-ä -dire en suivant la particule, s’ecrit : 

_ dVx _ V x {t + dt) - V x (t) 
üx dt dt 

La formulation d’Euler cherche ä exprimer l’acceleration en fonction des variations locales par rapport au temps et ä l’espace 
: pour distinguer les deux types de derivees, on notera ou encore une derivee locale par rapport ä t ou x. On n’utilisera 
la derivee droite que pour une derivee au sens clasique c’est-ä -dire Lagrangien. 

Relation entre acceleration Lagrangienne et Eulerienne Ecrivons dabord que V x (t + dt) = V x (t + dt) + dx . Ceci s’ecrit 
aussi: 

/ dV f)V 

V x (t + dt) = V x (t) + + ~Q^ dx 

d’oü l’expression de l’acceleration selon Ox : 


dVx _ &K dV^dx _ dV^ dV x 
dt dt dx dt dt x dx 


En 3 D on aurait ecrit : 


expression que l’on note aussi : 


- dVx t/ dVx , t/ dVx t/ dVx 

“ 777 1" Vx 1" Vy — h V z — — 

ot ox oy dz 


ßv -> - 
ax = ~gf + {V.V)V X 

On aurait trouve des expressions äquivalentes pour a y et a z de Sorte que la relation finale liant l’acceleration particulaire 
(Lagrangienne) et locale (Eulerienne) s’ecrit : 

et l’Equation d’Euler d’un fluide parfait soumis ä des forces de pression, ä la gravite et ä d’eventuelles forces externes fy 
s’ecrit : 


dV 

P ld =P 


dV_ 

~dt 


{V.S7)V 


= - gradP + pg + f v 


Le formulaire joint en fin de document donne les diverses expressions des Operateurs usuels en coordonnees cartesiennes 
ou cylindriques. 

Cas particuliers : 


• Ecoulement permanent : la derivee par rapport au temps vaut 0 , cad. que la vitessse en un point quelconque ne varie 
pas avec le temps soit ^ = 0. 

• Ecoulement etabli : signifie que l’ecoulement cesse de varier en espace. Lorsqu’un fluide penetre dans un tuyau par 
exemple, il faut une certaine longueur (de l’ordre de 30 ä 50 fois le diametre hydraulique) pour que le profil de vitessse 
ne varie plus suivant la direction d’avancement du fluide : on dit que l’ecoulement est alors etabli d’oü : yr- = 0. 


L’equation vectorielle d’Euler peut etre projetee sur les 3 directions de l’espace et fournit ainsi 3 equations. Toutefois, un 
Probleme Standard 3D de Mecanique des Fluides possede 4 inconnues : 3 de vitesse et 1 de pression : il manque donc une 
equation que nous allons voir maintenant. 
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2.4 Conservation de la matiere 


Un bilan de masse sur une surface S entourant un volume V doit traduire que, pendant un intervalle de temps dt, la difference 
entre les debits de masse entrant et sortant est necessairement egale ä la masse de matiere accumulee dans le volume V pen- 
dant le temps dt. Ceci s’ecrit: 


d_ 

dt. 


pdv + / p.(V .n e )da = 0 
i Js 


Le nrembre de gauche represente la Variation temporelle de masse i^jf ) du volume v: , tandis que le mernbre de droite 
represente le bilan global de masse ayant traverse la surface S entourant le volume v. 

Le theoreme d’Ostrogradski pennet de transformer cette relation : il pennet en effet d’evaluer une integrale de surface 
d’un champ de vecteur A ä partir de la divergence de ce champ de vecteur dans le volume enferme par la surface, soit: 


/ A.n dS= divA dv 

Js Jv 

oü n designe la normale orientee vers l'exterieur. Rappeions que la divergence est un Operateur qui s’applique ä un vecteur 
et traduit les variations du flux de ce vecteur au travers d’un volume infiniment petit place autour du point oü l'on calcule 
cette divergence. En coordonnees cartesiennes il s’ecrit: 


div A = 


dA r 


dAy dA z 


dx 


dy 


dz 


L’equation de Conservation de la matiere devient donc: 

[dp 


dt 


.dv + / div{pV)dv = 0 


qui, en faisant tendre le volume v vers 0 donne: 


^ + div (pV) = 0 


cette relation peut aussi s’ecrire sous une autre forme, äquivalente: 

dp 


dt 


+ p div V = 0 


Pour un fluide incompressible (p = Cste) en ecoulement permanent, cette relation devient plus simple et s’ecrit: 


divV = 0 
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Chapter 3 

Bilans de matiere et d’energie dans un 
ecoulement 


L’objectif de ce chapitre est d’apprendre ä faire un bilan ä travers une surface : deux types de bilan seront faits : des bilans 
sur des quantites scalaires (sans direction) telles que la masse, l’energie, un solute transporte par l’ecoulement ou encore 
des bilans de quantite de mouvement qui est une grandeur vectorielle : l’objectif ce travail etant d’exprimer le theoreme de 
quantite de mouvement dont l’aplication permet de quantifier la force qu’exerce un ecoulement sur une paroi. 


3.1 Bilan de matiere 


Cette expression est aussi appelee bilan de masse : on clierche ici ä chiffrer le debit massique ou volumique qui traverse une 
surface. 

Cas particulier : Champ de vitesse constant, de valeur V et perpendiculaire ä la surface S: 

Le debit volumique vaut simplement : Qv = VS et le debit massique Q m = pVS. 

Dans le cas le plus general (cf. Figure 3.1) les expressions precedentes deviennent alors : 


Qv 


V.ndS 


et 


Qn 


pV.ndS 



Figure 3.1: Debit d’un champ de vitesse ä travers une surface S 


- Si p = Cste alors : Q m = pQy - Bilan ä travers un tube de courant : cf. Figure 3.2 :on a alors dans le cas le plus 
general (p i ^ P 2 ) d’un ecoulement compressible : f f s p\V\.ndS = J f s piVi.ndS. Si p\ ~ P 2 et que V\ et V 2 ~ Cste alors 
: V 1 S 1 = V 2 S 2 = Qv- ~ Sinon, on peut toujours revenir ä la definition de la vitesse debitante V = ^ et ecrire : V\S\ = V 2 S 2 
Exemples de calcul de debit ä travers une surface: 

Exemple 1 : Debit d’un profil de vitesse parabolique en conduite cylindrique (cf Figure 3.3 ) 
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Figure 3.2: Debit ä travers un tube de courant 


r 



Figure 3.3: Exemple 1 : debit d’un profil parabolique en conduite cylindrique 


Le champ de vitesse est defini par : 


r 2 

V x (r) = V^ax(l - -j^) 

On en deduit le debit volumique: 


d’oü: 



V x 2n rdr = 2TrV max { — 


& 

T 


) 


v = 


v m 


Exemple 2 : debit d’un profil parabolique entre plans paralleles cylindrique (cf Figure 3.4 ) 



Figure 3.4: Exemple 2 : debit d’un profil parabolique entre plans paralleles 


Le champ de vitesse est ici defini par : 


V x (y) = V m 


■('-ff) 
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et on en deduit le debit volumique par unite de largeur : 


O r h h 4 

J= 2 J q V*dy = 2 V max (h --) = -hV m 


d’oü : 


— _ 2 
r — —Vmax 

Exemple 3 : Ecoulement tournant entre deux cylindres concentriques (cf Figure 3.5 ) : 



Figure 3.5: Exemple 3 : debit d’un ecoulement tournant entre deux cylindres 


Cette Situation se rencontre par exemple dans un viscosimetre, appareil servant ä mesurer la viscosite d’un fluide que Ton 
ne connait pas (c’est entre autres le cas d’un nouveau produit que l’on vient de formuler) : le cylindre interne tourne ä la 
vitesse angulaire oj tandis que le cylindre externe est ä l’arret. Si les rayons R\ et R 2 sont proclies, on peut montrer que le 
champ de vitesse s’exprime ainsi: 


Ve(r) 


ujRi 

R2 — Ri 


(i?2 


d’oii le debit tournant par unite de hauteur du cylindre: 


r) 


Q 

h 


d’oii la vitesse debitante : 


f R \ iJr= , Rl ^L 

/Hi Z 


- = ohRi 


qui est donc la vitesse moyenne entre la paroi interne qui tourne ä ujR± et la paroi externe qui est ä vitesse nulle 
Exemple 4 : Ecoulement axial entre deux cylindres concentriques (cf. Figure 3.6) : 

Si les rayons Ri et i ?2 sont tres differents on doit cliercher la solution du profil de vitesse V z (r) et l’integrer entre ces 
memes bornes soit : 

rR 2 

Q= / V z dr 
Jr 1 

Si i ?2 est tres proche de R\ de Sorte que l’onpeut ecrire : R- 2 ~ R\ + e alors onpeut se servil- du resultat de l’exemple 2 et 
on obtient : 


Q — ‘I'ItRg ^ Vmax 
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r 2 



Figure 3.6: Exemple 4 : debit d’un ecoulement axial entre deux cylindres 

3.2 Bilan de quantite de mouvement 

Intuitivement, tout le monde a deja ressenti la notion de ”force” d’un ecoulement oü d’un jet d’eau. On pourrait penser 
que cette ” force” est seulement liee au debit mais l’exemple de la lance ä incendie qui suit contredit ce sentiment. Dans la 
premiere experience, un jet d’eau ä 1 m/s est issu d’un tuyau horizontal et frappe un plan vertical rnonte sur un ressort. Le 
ressort se comprime de fagon ä equilibrer la force que le jet exerce sur la plaque. Dans la seconde experience, le tuyau est 
equipe d’un convergent de Sorte que la vitesse de sortie du fluide augmente. Pourtant le debit est le meme que dans le cas 
precedent. On voit alors que le ressort est plus comprime, traduisant le fait que la force exercee par l’ecoulement sur la plaque 
est plus elevee. On vient donc de montrer que, en plus d’un debit massique ou volumique, un ecoulement possede aussi un 
certain debit de quantite de mouvement. Un autre exemple illustrant la force d’un ecoulement est la technique industrielle 
de decoupe de materiaux metalliques qui consiste ä utiliser un jet d’eau sous tres haute pression (plusieurs centaines de bars). 




Figure 3.7: Illustration de la quantite de mouvement d’un ecoulement. 


La quantite de mouvement par unite de volume est le vecteur pV (en N.m~ 3 ). Comme beaucoup d’autres quantites 
scalaires ou vectorielles, la quantite de mouvement est une grandeur qui est transportee par l’ecoulement et pour laquelle on 
peut faire un bilan. Par analogie avec ce qui a ete montre pour calculer un debit de masse, le debit de quantite de mouvement 
(i qdm ) ä travers une surface £ s’ecrit: 


pV(V.n)da 


Considerons maintenant une surface £ fermee, entourant un volume v : le theoreme de quantite de mouvement exprime 
sur cette surface s’ecrit: 
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F e = J pV{V.fi)da 

• Le terme de gauche represente la contribution des forces externes qui peuvent diminuer ou augmenter la qdm du fluide 
contenu dans le volume v. Ce sont: 

• les forces de pression (surfaciques) qui peuvent diminuer ou augmenter la qdm 

• les forces de viscosite qui, le plus souvent, consomment de la qdm 

• la force de pesanteur 

• les forces de volume qui peuvent etre electriques, thermiques etc... 

• Le terme de droite, cornrne on l’a vu, represente le bilan global au travers de la surface fermee E qui est souvent 
denommee "Surface de controle”. 

Protocole d’application du theoreme de quantite de mouvement: 

- choisir un Systeme d’axes Orientes 

- choisir une surface fermee adaptee 

- orienter les normales ä cette surface vers l’exterieur 

- projeter les equations de qdm sur chaque axe 

Exemple 1: Ecoulement etabli en conduite cylindrique (cf Figure 3.8) 



D 

< 

Figure 3.8: Calcul de la contrainte sur un tube. 



L’ecoulement etant etabli signifie que le champ de vitesse ne varie plus selon Ox, d’oü = 0. En choisissant la surface 
ABCD (qui est un tronc de cylindre) on en deduit J f Sl+ g 2 pV x (V.n)dS = 0. D’oü : +P\ttR 2 — P 2 TtR 2 — tp2ttRL = 0 ce 
qui donne la valeur de la contrainte parietale: 


Tp 


(■ Pi-Pi)R Al -2 

2 L 4 2 pV 


Exemple 2 : Poussee sur un coude (cf Figure 3.9) 

Valeurs numeriques : P\ = 16.10 5 Pa, S\ = 0.08?n 2 , V\ = 5m/s , : P 2 = 10.10 5 Pa, S 2 = 0.01?7? 2 , pf = 1000 kg/m 3 
Posons : R: action du fluide sur le coude d’oü: — R: action du coude sur le fluide. En negligeant le poids, la projection selon 
Ox et oy donne: 


-pVfSx = PiSi - R x 
+ P ViS 2 = -P2S2 - Ry 


d’oü: 


R x = P 1 S 1 + pV/S-i = P 1 S 1 + pQV 1 = 130000V 
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c 


D 



R y = -P 2 S 2 - pViS 2 = -P 2 S 2 - pQV 2 = -26007V 

d’oü l’angle que fait la force avec l’horizontale:# = Arctg(R y /R x ) = 11, 3° 

Exemple 3 : Poussee d’un jet axisymetrique liquide sur une plaque (cf Figure 3.10) 



! U 


Figure 3.10: Calcul de la poussee sur un coude. 


Hypotheses : la pression P est contante dans tout le jet liquide (P = P a tm ) On constitue une surface fermee par 

S e + Sl + Sint 

On appelle R l’action du jet sur la plaque. La projection sur Ox donne : 

-pV 2 S = P 0 S int - R x 

Bilan sur la face interne de la plaque : R x = pV 2 S + PoSi n t 
Bilan sur la face externe : — PoSi n t 
D’oü la resultante : +pV 2 S. 

Exemple 4 (ä resoudre) ijonction de 2 ecoulements (cf Figure 3.11) 


Dans la section A 2 ecoulements de vitesse et de section distinctes se melangent : 

- Determiner V 2 = f{V\) 

- Exprimer Pa — Pb en fonction de (p, V \ ) 
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A B 



Figure 3.11: Jonction de deux ecoulements. 
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Chapter 4 

Dynamique des fluides visqueux 


Nous avons vu que l’equation d’un fluide parfait (non visqueux) soumis ä une force volumique f v incluant la gravite est 
donnee par l’equation d’Euler: 

dV -w 

P- 7 T = -grad p + } v 
dt 

Comme on peut le pressentir, la viscosite cree dans un ecoulement des contraintes, denonnnees contraintes visqueuses, qu’il 
convient de prendre en compte si Fon veut etablir l’equation d’un fluide reel visqueux : l’objet de ce chapitre est d’etablir 
l’expression du terme traduisant ces effets visqueux et qu’il faudra donc ajouter ä l’equation d’Euler pour la completer. 
Avant cela donnons quelques exemples qui illustrent l’effet de la viscosite. 

Elargissement d’un jet :cet exemple se rencontre par exemple en combustion lorsqu’on injecte sous forme de jet le carbu- 
rant dans la chambre de combustion. Dans le cas des moteurs d’Ariane le melange 0 2 — H 2 se fait par deux jets cylindriques 
concentriques (cf. Figure 4.1). On observe l’effet de lissage de la viscosite qui tend ä homogeneiser le champ de vitesse qui 
de ce fait s’elargit rnais diminue d’intensite. On retient donc que la viscosite tend ä gonnner les differences de vitesse. 

Le second exemple (cf. Figure 4.2) illustre la celebre experience de Reynolds dans laquelle on injecte un filet ebenere dans 
un fluide au repos. Pour une tres faible vitesse d’injection, c.a.d. un faible Reynolds, l’encre se fraie un clremin en diffusant 
tres peu dans la direction transverse. Lorsqu’on augmente la vitesse d’injection, la dispersion devient alors majeure. 


4.1 Experience de Couette 

Cette experience (cf. Figure 4.3) sert ä quantifier l’effet de la viscosite lorsqu’on varie le champ de vitesse. Plus precisement 
c’est surtout l’influence du cisaillement que Fon cherche ä quantifier c.a.d. l’influence des variations spatiales de vitesse 
(gradients). Le fluide etliche remplit l’espace situe entre deux cylindres concentriques. Le cylindre interne tourne ä la vitesse 
uj tandis que le cylindre externe est fixe. La mesure du couple mecanique resistant sur l’arbre moteur donne acces ä la 
valeur de la force tangentielle F qui freine le cylindre interne. Appelons V la vitesse lineaire de la paroi interne, S la surface 
mouillee, e l’entrefer, et F la force resistante au mouvement: 

L’experience montre que le rapport: 


F 

sTC 


est constant pour un grand nombre de fluides usuels (air, eau,...). Cette constante est la viscosite dynamique fi. On voit 
donc que la contrainte r, rapport de la force F ä la surface mouillee S, s’ecrit: 


T 


/ 

s 



avec: V = ojR et S = 2t: R h. La quantite V/e represente en fait le gradient de vitesse ^ oü r represente la direction 
normale ä F ecoulement. 

Pour cette experience particuliere, on voit donc que la contrainte r vaut: r = Cette contrainte s’exprime en 

N/m 2 . Pour un ecoulement tridimensionnel quelconque d’un fluide incompressible, il existe en fait un champ de contraintes 
beaucoup plus complexe que celui de l’ecoulement de Couette. II faut donc definir d’abord un Systeme de notation rigoureux. 
Considerons la figure 4.5 oü r,;j designe la contrainte visqueuse de direction i exercee sur une face de normale j. 

On montre qu’un grand nombre de fluides voient leur contrainte definie par: 
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Figure 4.1: Dispersion du champ de vitesse dans la chambre de combustion d’un moteur de fusee 




Fluide visqueux 


Encre 
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\ 
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Figure 4.2: Experience de Reynolds 
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Figure 4.3: Experience de Couette 


Tij = ß 


( dUi 


dUj \ 

dxi J 


Cette loi de comportement est celle des Fluides Newtoniens . 

Toujours pour un fluide incompressible, on tire facilement les expressions des contraintes en fonction des composantes 
(u, v , w ) de la vitesse : 


du 

T~xx — 

dv 

T yy = 


= 2/i 


dy 

dw 

dz 




( du 


dv' 

T yx - 

= y 1 

\dy 

+ 

dx 



f du 


dw 

T zx ~ 

= M( 

ydz 

+ 

dx 


1 

f dv 


diu 

T zy = 

= ld \ 

K dz 

+ 

dy 


Lorsqu’on cherche ä faire le bilan des contraintes sur un volume de contröle, il convient d’ exprimer la contrainte en 
fonction de la normale n ä chaque face. On ecrira: 

dV 

T = ^ 

ou n designe la normale dirigee vers l’exterieur. Considerons le volume de contröle ABCD de la Figure 4.4 inclus dans un 
ecoulement unidirectionnel variant dans la direction perpendiculaire ä BC et AD. La contrainte qui s’exerce sur BC aura 
un signe > 0 du fait que ^ > 0 (la vitesse augmente quand on se deplace dans la direction ni). La contrainte s’exercant 
sur AD sera < 0 car ^ > 0 (la vitesse diminue quand on se deplace dans la direction n. 2 ). 


Ainsi, si le profil de vitesse est lineaire, la quantite ^ sera egale, en valeur absolue, sur BC et sur AD, de Sorte que la 
force resultante sur le volume ABCD sera nulle. 


4.2 Equations de Navier- Stokes 

Nous venons de voir comment calculer le champ de contrainte r en fonction des variations spatiales (des gradients) du champ 
de vitesse. Toujours dans l’objectif de determiner le terme du ä la viscosite que nous souhaitons ajouter ä l’equation d’un 
fluide parfait (equation d’Euler) il faut maintenant determiner la force volumique qu’exerce un champ t de contrainte existant 
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dans le fluide. 


Considerons le volume elementaire represente sur la figure 4.5. 

De meine que la pression n’intervenait que par son gradient, on montre que les contraintes visqueuses n’interviennent 
que par leur divergence et produisent, par exemple suivant Ox le terrne suivant : 

3t 3t 

w 1 xx . w ' xy . w ‘ xz 

dx dy dz 

C’est donc par sa divergence que le tenseur des contraintes visqueuses r interviendra dans l’equation du mouvement d’un 
fluide reel. 

Par unite de volume, et toujours dans la direction x, l’equation du mouvement d’un fluide reel s’ecrit donc: 


dU x 

P ~dT 


dp 

dx 


+ div t x + f x 


oü f x designe les forces de volume s’exercant sur le fluide dans la direction x. 
Une ecriture vectorielle plus compacte de cette equation sera: 


av -> _ -* 

p—rr = -Vp + divr + f v 
dt 

En reportant les relations ci-dessus dans l’equation du mouvement, on obtient, pour la composante x: 

du dp d 2 u ( d 2 u d 2 v \ f d 2 u d 2 w \ 

P dt dx P dx 2 P \ dy 2 dxdy ) ^ \ dz 2 ^ dxdz ) + x 

Apres avoir remarque que la somme de trois des termes de cette equation est nulle ( car egale ä d(divV)/dx), et apres 
avoir mene des calculs similaires pour les deux autres composantes y et z, on obtient l’equation vectorielle suivante: 

dV 

P~TT — ^ Vp + pV V + fv 

dt 

C’est l’equation de Navier-Stokes. Dans cette formulation, f v s’exprime en N/m 3 et p = pv represente la viscosite 
dynamique du fluide. Le formulaire place en Annexe fournit les equations de Navier-Stokes en coordonnees cartesiennes et 
cylindriques. 


4.3 Les conditio ns aux limites 

Les composantes de l’equation de Navier-Stokes fournissent un ensemble d’equations differentielles auquel il convient d’adjoindre 
les conditions aux limites adequates. 

• Pour une paroi: la condition d’adherence impose: 

V = 0 




Pour un axe de symetrie 

En cartesien: Prenons l’exemple de l’ecoulement entre deux plans fixes paralleles ä Ox. 
ecrira: 


dv x 

dy 


= 0 et v y = 0 


Sur l’axe Ox de symetrie on 


*En coordonnees cylindriques et pour un ecoulement axisymetrique, on ecrira, sur Taxe: 


u r 


3 U 

0 et ug = 0 et — - ■ 0 
dr 


• Pour une surface libre: 


n + dUx 
Uy = 0 et — — 

dy 


du z 

dy 


= 0 
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Figure 4.6: Differentes conditions de symetrie 


4.4 Viscosite des fluides non-Newtoniens 

Beaucoup de fluides s’apparentent ä un fluide Newtonien. Toutefois, il n’est pas toujours possible de les considerer comme tels. 

Nous avons vu qu’un fluide Newtonien est tel que la contrainte r qu’il exerce est proportionnelle au taux de deformation. 
En ecoulement unidirectionnel, cela s’ecrit: 

dU x 

T= »w 

oü h est une constante caracteristique du fluide. 

Pour un fluide non-Newtonien, on pose quand-meme: 


T 


dU x 

dy 


(4.1) 


oü fi a est la viscosite apparente du fluide qui ne sera donc plus une constante. On distingue plusieurs classes de fluides 
tels que: 


• Fluides de Bingham 

La relation contrainte-taux de deformation s’ecrit: 


dU x 

T = To + Ho —j 

dy 

On note dans ce cas l’existence d’un seuil de contrainte au dessous duquel le fluide ne se met pas en mouvement (cf 
courbe numero 4). 

• Fluides d’ Ostwald 

Leur loi de comportement se met sous la forme: 

r = C(du/dy) n = C(du/dy ) n ~ 1 du/dy 
par identification avec la formule 4.1 on obtient: 


Ha = C(du/dy) n 1 


Si n = 1, le fluide est Newtonien (courbe 1) 

Si n > 1, le fluide est dit dilatant: sa viscosite augmente avec le taux de deformation. On dit aussi qu’il est rheo- 
epaississant. Exemples: amidon, guimauve (courbe 5). 

Si n < 1, le fluide a un comportement pseudo-plastique, c’est-ä-dire que sa viscosite diminue avec le taux de deformation 
(margarine, huile vegetale, mayonnaise, ...). On dit que le fluide est rheo-fluidisant (courbe 6). 


27 


dy 

Fig. 1 5. - Exemples de lois de comportement 

1 . Fluide newtonien - 2. Fluide parfait - 3. Solide parfait - 4. Fluide de Binghair» - 5. Fluide 
dilatant - 6. Fluide pseudo-plastique. 

Figure 4.7: Exemples de lois de comportement d’apres Comolet 
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Chapter 5 


Ecoulements en conduites et canaux 


Les exemples qui suivent illustrent des situations d’ecoulements laminaires simples en conduites pour lesquels une solution 
analytique peut facilement etre obtenue. 

5.1 Ecoulement d’un fluide visqueux incompressible entre deux plans par- 
alleles 

Considerons l’ecoulement d’un fluide visqueux s’ecoulant entre deux plaques paralleles infinies sous fleffet d’un gradient de 
pression constant. 



Figure 5.1: Ecoulement de Couette entre plans paralleles 


Hypotheses: Supposons l’ecoulement permanent ( = 0 ), plan = 0 ) et infini dans la direction x = 0 ). La gravite 
est dirigee suivant Oz. La geometrie inscite ä cherclier une solution de la vitesse du type (m(ä), 0, 0). L’incompressibilite 
irnpose: 


du 



tandis que les trois composantes de NS s’ecrivent: 

1 dp d 2 u 

suivant x : — — = v-^—r 
p dx dz 2 

• f 1 d P n 
suivant y : — — = U 
pdy 

■ r 1 dp n 

suivant z: — g 

p dz 

d’oü, en introduisant la pression etoilee p* = p + pgz: 


dp* 

dz 


On en deduit l’equation differentielle suivante pour u: 


0 


d 2 u 1 dp 

dz 2 p, dx 
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d’oti: 


Les conditions aux limites: 


1 dp z 2 

u = -- — - + Az- 
pdx 2 


B 


u = 0 en z = ±- 


h 


fournissent la solution finale: 


l dp f 4 z 2 — h 2 


Le profil de vitesse est donc parabolique. On notera que, pour que l’ecoulement soit dirige vers les x Croissants, il faut 
que le gradient de pression soit < 0 ( la pression decroit suivant x). On peut noter les points suivants: 


la vitesse maximum en z = 0 est teile que: 


h 2 dp 

^ max ö ö 

8/i OX 

• le debit volumique par unite de largeur des plaques vaut: 


Q s’exprinre en m 3 /s/m. 

• la vitesse moyenne V vaut: 


C+h/2 



u{z)dz = 


h 3 dp 
12/r dx 


• la contrainte t xz aux parois vaut: 


y = Q = tf_d]>_ = % u 

h 12 pdx 3 U 


du dp h 
T “ = ~ ß d~z = *^'2 “ 


z=+h/2 et -h/2 


Le signe - dans l’expression de la contrainte vient du fait que la contrainte sur la paroi est egale et opposee ä celle subie 
par le fluide. 

* Cas oü la plaque superieure est animee d’une vitesse uniforme: dans ce cas, le glissement de la plaque (qu’on supposera 
de vitesse V dirigee vers les x > 0) peut etre eventuellement conjugue avec un gradient de pression externe. Appelons G ce 
gradient de pression (G = ff )• Si G < 0, le gradient de pression est moteur et amplifie le phenomene cree par la plaque. Si 
G > 0 le gradient de pression est dit adverse et s’oppose ä l’action de la plaque. 

i z 

| vitesse V 


I 



Figure 5.2: Ecoulement de Couette avec une paroi mobile 


Les equations sont les memes que dans le cas precedant, seules les conditions aux limites changent. Elles s’ecrivent: 


On obtient alors la solution: 


u = 0 

en z = —h / 2 

u = V 

en z = +h/2 


G , A 2 Vz V 
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• le debit par unite de largeur vaut maintenant: 


f+h/2 y h 

Q= u(z)dz = — 

J-h/2 1 

• Si G = 0 u = yZ + tt , et la repartition de vitesse est lineaire. 

• la contrainte tangentielle vaut ici: 


G 

12 Jl 


h 3 


Txz = - 


du\ 

[ " dz )hl2 


avec 


du 

dz 


Gz 


V 

h 


soit: 




V . Gh 

Y 


5.2 Ecoulement d’un fluide visqueux dans une conduite cylindrique 


Soit une conduite cylindrique de rayon R dans laquelle s’ecoule un fluide visqueux sous l’action d’un gradient de pression 
G = ^ (on utilise maintenant les coordonnees cylindriques r , 9, z. 

Hypotheses on suppose l’ecoulement axisymetrique = 0) et permanent(Jj = 0 ). Cherchons une solution du type 
{u r = 0; ug = 0; u z = u{r)). 

Remarque la gravite n’est pas prise en cornpte dans ce calcul: on sait, qu’etant perpendiculaire ä l’ecoulement, eile ne 
fait que modifier le champ de pression. 



Z 


Figure 5.3: Ecoulement de Poiseuille en conduite cylindrique 


La continuite impose: 

tandis que les trois composantes de NS s’ecrivent: 


dz 


dp 

suivant r : — — = 0 
dr 


dp 

suivant 9 : — = 0 

d9 


suivant z : 


dp 

dz 


1 d ( du 


dr V dr 


= 0 


Les conditions aux limites « = 0 en r = Ä et §^=0enr = 0 donnent la solution suivante: 

ör 

u — — (r 2 — R 2 ) avec G < 0 
4/x v ' 
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• Calcul du debit: 


soit: 



R 2 ) 2i xrdr 


Ceci est la formule de Poiseuille. 

• Vitesse maximuni: en r = 0 on trouve: 

• Vitesse moyenne (ou encore vitesse debitante): 




Urn n er 


GR 2 


u = 


Q 

nR 2 


GR 2 

8 n 


Umax 

2 


Valeur de la contrainte: seule la contrainte T rz n’est pas nulle et vaut: 


du 

T rz = T zr = PVT~ = 

ar 

Notemment en r = R et sur la paroi, on trouve: 

GR 



Gr 


Application au viscosimetre: en faisant circuler dans une conduite circulaire de petit diametre un fluide visqueux on 
peut, ä partir de la connaissance du gradient de pression et du debit, en deduire la viscosite du fluide. Cela constitue une 
application de la loi de Poiseuille. 

Coefficient de perte de Charge en conduite cylindrique 


V 


V 


- i h 




Diametre D Q 

Longueur L 

Figure 5.4: Mesure de la perte de Charge en conduite cylindrique 


En negligeant les pertes de Charge singulieres, appliquons le tlieoreme de Bernoulli entre l’entree et la sortie du tuyau: 

L V 2 L V 2 

H t = H‘2 + X > X = h 

D 2g D 2g 

En remplacant l’expression de V par sa valeur en fonction de G et G par: G = on deduit l’expression du coefficient 

de perte de Charge A: 

64 
~ ~Re 

oü Re designe le nombre de Reynolds base sur le diametre de la conduite: 

Re = YR 
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5.3 Ecoulement de Couette cylindrique 

Considerons le cas oü seule la paroi externe tourne. La encore nous nous limiterons ä un regime stationnaire et axisymetrique. 
Clierchons un champ de vitesse de la forme: (0; v(r); 0). 



Figure 5.5: Ecoulement visqueux entre deux cylindres tournants 


La composante radiale de l’equation de Navier-Stokes traduit l’equilibre entre le gradient radial de pression et la force 
centrifuge: 

1 dp v 2 

p dr r 

La composante suivant 0 s’ecrit: 

v 

Vw - ~ä = 0 

Suivant z, on verifie l’equilibre hydrostatique: 

P*(r) =P~ pgz 

C’est donc bien la resolution de l’equation azimutale qui va donner la solution de l’ecoulement. Pour cela, clierchons une 
solution de la forme r n . En introduisant ce type de solution dans l’equation de v on trouve: 

„n— 1 y,n 

n{n — l)r” -2 + n j = 0 


d’oii l’equation de n: 
dont les Solutions sont: 


n(n — 1) + n — 1 = 0 


n = ±1 


la solution generale de v s’ecrit donc: 


B 


v(r) = Ar H 

r 


Conditions aux limites : v(R\) = 0 et u(i? 2 ) = 

Ceci fournit la solution finale: 

VtR 2 2 __ VtR\Rl 1 
V ~ Rl - R\ ' r ~ R 2 — R\ r 

• Champ de pression: la pression etoilee p* (r) s’obtient par Integration de pv 2 

• Couple resistant: la contrainte tangentielle a r g s’ecrit: 

( dvg Vg\ 

<?r8 = ß\ soit en R '2 : (T r g = - 

\ dr r J 

On en deduit le couple C de freinage de la paroi externe: 


/r 
2 pB 

'“rT 


d2d2 

C = F.R 2 = a r g2TrR%h = 

R‘2 — Ri 

• Cas des petits entrefers: En posant Ri = R et R 2 = R + 6 avec 5 « R, on peut faire un developpement limite de v 
qui, naturellcment, devient lineaire entre les deux parois. 
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Chapter 6 


Ecoulements sur paroi et autour d’obstacles 


Par Opposition aux ecoulements ä l’interieur de conduites, encore appeles ecoulements internes, ce chapitre traite le cas 
d’ecoulements dits ” externes ” , s’ecoulant autour d’obstacles. On considerera d’abord le cas de la plaque plane qui pennet 
de definir la notion de couche limite, puis on considerera le cas d’ecoulements autour de cylindres, spheres et autres geometries. 

6.1 Ecoulement sur plaque plane 

6.1.1 Notion de couche limite 

Dans beaucoup de situations, et notamment dans le cas d’ecoulements au voisinage de parois solides, la viscosite joue un role 
fondamental qu’il convient d’examiner avec plus de detail. Au voisinage d’un obstacle, par exemple, la condition d’adlierence 
ä la paroi impose u = 0 ä cet endroit. Au loin, par contre, F ecoulement est anime d’une vitesse externe, ce qui va donc 
donner naissance ä une zone de raccordement, dite Couche Limite dans laquellc les gradients de vitesse sont eleves. Dans 
cette region, F ecoulement est dit cisaille . Comrne nous allons le voir dans le paragraphe qui suit, on rencontre aussi des 
ecoulements cisailles dans des situations oü il n’y a pas de parois. 

Exemples de couches cisaillees dans un fluide visqueux : 

• Couche limite au voisinage d’une paroi 

Comrne on le voit sur les exemples de la figure 6.1, l’existence d’une paroi se traduit par un inflechissement du profil 
de vitesse qui va jusqu’ä s’annuler sur la paroi elle-meme. L’epaisseur de la couche limite va croissante au für et ä 
mesure de la progression du fluide le long de l’obstacle. Le fluide perd de la quantite de mouvement dans cette zone 
sous l’effet des contraintes visqueuses qui sont beaucoup plus elevees ä cet endroit que dans le reste de F ecoulement 
dit ” ecoulement externe”. Aussi bien d’un point de vue descriptif que mathematique, on est amene ä distinguer deux 
zones distinctes: 

— la couche limite proprement dite oü les gradients de vitesse longitudinaux sont beaucoup plus faibles que les 
gradients de vitesse transversaux. Cette hypothese, dite approximation de la couche limite, est a l’origine de 
grandes simplifications dans le traitement mathematique des equations de Navier-Stokes dans cette zone. 

— F ecoulement externe oü la viscosite intervient tres peu et oü l’ecoulement peut etre condiclere comrne celui d’un 
fluide parfait. L’une des contributions essentielles de l’ecoulement externe est d’imposer son charnp de pression 
ä la couche limite. Ainsi, on peut supposer que le gradient de pression transverse dans la couche limite est tres 
faible tandis que le gradient longitudinal est impose par l’ecoulement externe. 

• La couche de melange 

Lorsque deux ecoulements de vitesses differentes se rejoignent, comrne le rnontre la figure 6.2, il y a formation d’une zone 
cisaillee oü s’effectue un lissage progressif de la discontinuite de vitesse. Cette zone de faible epaisseur, dite "Couche 
de melange” donne egalement naissance ä de grands gradients transverses de vitesse. Son epaisseur augmente lorsque 
l’on se deplace vers l’aval, sous l'effet de diffusion transverse qu’assure la viscosite. 

• Les jets 

le jet est un ecoulement cisaille cree par l’expulsion d’un fluide dote d’une certaine quantite de mouvement dans un 
fluide au repos. Comrne on le voit sur la Figure 6.3-a, une zone de melange se forme ä la frontiere du jet et, plus en 
aval, ces deux zones se rejoignent pour former un profil etabli. Ce jet peut eventuellement etre confine (Figure 6.3-b) 
comrne dans les injecteurs coaxiaux par exemple. 
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Couche limite au voisinage d'un profil plongi 
dans un ecoulement. Pour plus de clarti, 
l’epaisseur de la couche limite est exagirie. 
Dans la plupart des cas pratiques, l'ipaisseur 
de la couche limite est de quelques mm. 



Couche limite sur une plaque plane. 



Couche limite autour d'un 
corps de forme conique. 




Couches limites formies par l'äcoulement initial dans un conduit Apris une 
certaine distance, les de ux couches limites se rejoignent et l'ecoulement 
dans le conduit est pleinement developpi. 


Figure 6.1: Quelques types de couches limites de paroi. 



Figure 6.2: Couche de melange ä la confluence de deux ecoulements parallelles de vitesses differentes. 
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Figure 3a). Jet libre formi en sortie 
d'une tuyire. 



Figure 3b). Jet confini ä l'intirieur 
d'un ijecteur. 


Figure 6.3: Formation d’un jet libre a)libre b)confine. 


• Les sillages 

Lecoulement derriere un obstacle engendre un deficit de vitesse comme le nrontre la figure 6.4. Cette discontinuite a 
eile aussi tendance ä se lisser lorsque l’on s’eloigne de l’obstacle, toujours sous l’effet de la diffusion transverse assuree 
par la viscosite. L’ensenrble de la zone perturbee constitue le sillage. 


U 



Figure 6.4: Sillage derriere un obstacle: ä l’interieur du sillage, la vitesse passe par un maximum. 


6.1.2 Coefficient de frottement sur plaque plane 

Des que l’ecoulement aborde la plaque, la couche limite se construit et s’epaissit au für et ä mesure que l’abscisse x mesuree 
le long de la plaque augnrente. L’epaisseur <5 de la couche est donc une fonction croissante de x. Ainsi, l’epaisseur qui separe 
la paroi de l’ecoulement dit ” externe ” oü la vitesse est globalement egale ä celle du fluide arnont, augnrente. De ce fait, 
la contrainte de frottement diminue avec x, en partant d’une valeur theoriquement infinie en x = 0 puisque le gradient de 
vitesse tend lui aussi vers l’infini. On construit un coefficient de frottement sans dimension, rapport de la contrainte parietale 
reelle Tp ä une echellc de contrainte valant \pU 2 : 



On trouve ainsi (calcul hors programnre) le coefficient de trainee adinrensionnel local (fonction de a:) : 

Cf{x) = 0 M4.Re~ h 
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oü Re x designe le nombre de Reynolds local: Re x = 

Coefficient de trainee: De la contrainte Tp, on deduit, la trainee D (D pour Drag en anglais, exprimee en Newton) qui, 
pour une plaque de largeur b et de longueur L s’ecrit: 


r 1 i 

D = I C f-pU 2 bdx 

Jo 2 


On rend cette valeur sans dimension en introduisant un coefficient de trainee Cp defini cornrne le rapport de la trainee 
D ä (l/2)pU 2 S oü S = bL designe la surface de la plaque. On obtient: 


D f‘ Cf yu 2 bdx 0.664 f‘ -Jx 
D \pU 2 S \ P U 2 bL L 

soit: 

C D = 1.328 Re^ 

Cette loi n’est valable que dans la zone laminaire. Lorsque l’abscisse x augmente, l’ecoulenrent dans la couche limite 
devient instable puis turbulent (cf. Figure 6.5). La valeur du nombre de Reynolds, construit avec la longueur x connne 
longueur caracteristique, qui separe le reginre stable du reginre instable est appelee valeur critique. Pour les ecoulements sur 
plaque plane cette valeur vaut: 


(Re 


. Ux 
Cc = — 


5.10° 


A la transition entre la zone laminaire et la zone turbulente, on observe une augmentation brutale de l’epaisseur de la 
couche limite (figure 6.6) ainsi que de la contrainte parietale. 



Figure 6.5: Transition entre le reguime laminaire et turbulent dans une couche limite sur plaque plane (d’apres Solliec, Cours 
de l’Ecole des Mines de Nantes). 


Pour un ecoulement turbulent sur plaque plane, le coefficient de frottement local s’ecrit (Figure 6.7): 

C f (x) = 0.0594i?e“ 1/5 

oü Cf est defini, cornrne dans la zone laminaire par: 

rp 


Cf = T 


PU 2 


On en deduit le coefficient de trainee Cp dans cette zone: 

Cd = ° 


\ P U 2 S 


= 0.074 Re 


- 1/5 


Cette formule est valable pour 5.10 5 < Re x < 10 7 . Pour Re > 10 7 on prendra plutot: 

_ 0.455 

c f = 


(log w Re x ) 2 - 58 

qui collc mieux aux resultats experimentaux et qu'il faudra integrer pour deduire la trainee totale D. 
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10 3 Cf 

Fig 7. Coefficient de 
frottcment de la plaque 
plane 

-..Cf = 0,05941?J 1/5 
Cf = 0, 0368 Rx 1/6 

Resultats exptrimentaux 
(ngmupts par Michel) 

A Dhauian 

1 O Schuttz-Grunow 

• Kcmpf. 



Figure 6.7: Coefficient de frottement local Cf(x) en fonction du nombre de Reynolds local Re x pour un ecoulement turbulent 
sans gradient de pression sur plaque plane. 
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6.2 Ecoulement autour d’obstacles 

6.2.1 Trainee et portance 

Considerons un ecoulement suivant un axe Ox (suppose horizontal) s’ecoulant autour d’un corps. Appelons Oz la direction 
normale (verticale dans ce cas). La force de resistance suivant Taxe Ox a dejä ete discutee et s ’appelle la trainee. Designons 
par F x la force correspondante. 

Si l’obstacle n’est pas symetrique par rapport ä l’axe Ox, il peut s’en suivre un effort non nul suivant Taxe Oz : designons 
par F z cet effort denomme portance. 

Trainee et portance sont les parametres des de tout objet volant, et meine dans une certaine limite, des parametres 
importants de dimensionnements de certaines voitures de courses (une portance excessive amene le vehicule ä s’envoler !). 
Les coefficients de trainee et de portance sont souvent exprimes sous forme adimensionnlle : 

c — ^ x 

x ~ \ P U 2 S 

C =-^_ 

* \pU 2 S 

oü S designe la surface frontale de l’obstacle (projetee suivant Ox). 

6.2.2 Decollement des couches limites 

Lorsqu’un fluide s’ecoulant suivant Ox le long d’une paroi est soumis ä un gradient de pression ’Tnoteur” c’est-ä-dire tel que 
dp/dx < 0 le fluide n’est pas ralenti et la couclie limite n’a pas de difficulte ä adlierer ä la paroi de l’obstacle qu’elle contourne. 
Ceci est illustre sur la figure 6.8 qui montre un ensemble” convergent-divergent” : dans la partie amont, la contraction de 
section entraine une diminution de pression(Bernoulli) et genere donc un grad p < 0. 



Figure 6.8: Ecoulements decolles dans des conduits fortement divergents: a) decollement symetrique b) asymetrique. 


En amont, 1’ ecoulement reste au contact de la paroi. Passe le col, la section augmente et le gradient de pression devient 
adverse ce qui entraine le decollement de la couclie limite et la creation d’une zone de recirculation. La figure 6.9 montre 
une vue grossie de la zone de decollement et l’on y voit que le point de decollement est caracterise par: 

du ( du\ . 

— = 0 — > t w = /U — =0 (contramte =0) 

9y \dyJ y =o 

6.2.3 Ecoulement autour d’un cylindre pour divers Re. 

Le decollement d’un ecoulement s’observe aussi autour d’un cylindre oü, au voisinage de 0 sa 90 degres l’ecoulement venant 
de l’amont decolle pour donner naissance ä deux tourbillons ”attaches” ä l’arriere du cylindre (figure 6.10). 

Ceci a lieu tant que le nombre de Reynolds, base sur le diametre du cylindre, reste inferieur ä Re ~ 40. Au dela, la 
couche limite devient instable et ses oscillations donnent naissance ä des tourbillons laches alternativement de chaque cote 
du cylindre. A l’aval, on observe une "Allee de tourbillons” encore denommee "Allee de Karman” (Figure 6.11). Quand 
Re augmente encore, le sillage devient turbulent et le point de decollement est deplace vers l’aval. Le nombre de Strouhal 
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Figur * il. Au point de ddcollement 
le profil de vitesse präsente une 
targente verticale et la contrainte 
pditale s'annule. 


Figure 6.9: Illustration du decollement d’une couche limite. 



Re > 40: Rue de tourbillons de Von Karman 

point de Separation 

P 0 “ 1 * ^ de la couche limite 



Re <3 «10* Re > 3 * 10 5 


Figure 6.10: Ecoulement autour d’un cylindre. 





S = tjj- caracterise la frequence adimensionnclle de lache des tourbillons. Les Figures 6.12 et 6.13 montrent Revolution du 
nombre de Strouhal S et du coefficient de trainee Cd avec Re. 



Figure 6.11: Allee de Karman generee par l’ecoulement autour d’un cylindre ä Re = 140 (d’apres Van Dyke: An Album of 
Fluid motion). 
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Figure 6.12: Nombre de Strouhal en fonction de Re pour V ecoulement autour d’un cylindre. 


6.2.4 Ecoulement autour d’obstacles quelconques 

La Figure 6.14 donne les coefficients de frottement d’ecoulements autour d’obstacles divers. 
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Figure 6.13: Coefficient de trainee de l’ecoulement autour d’un cylindre. 
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Figure 6.14: Coefficients de trainee pour un ecoulement autour d’obstacles divers - (D’apres Comolet) Cd = 


D 

±pU 2 S 
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Chapter 7 

Exercices 


Hydrostatique 

* (Schaum 2.14 p.21): Calculer la pression manometrique en A ( en bar) due ä la denivellation du mercure, de densite 
13,57 dans le manometre en U represente sur la Figure 7.1: 



Figure 7.1: 


Reponse : Pa = 1.0058 bar. 


* (Schaum 2.26 p.29):Trouver la difference de pression entre A et B dans le montage de la Figure 7.2: 
Reponse : Pa — Pb = 13.3 kPa. 


* (Schaum 2.41 p.34) Se reportant ä la Figure 7.3 pour quelle valeur marquee au manometre A la glycerine s’eleverat-elle 
au niveau B ? Les rnasses volumiques de l'huile et de la glycerine sont respectivment de 832 et de 1250fcg/?77. 3 . 

Reponse : 0.356ar. 


* (Schaum 3.3 p. 38) Calculer la force resultante P due ä l’action de l’eau sur la surface rectangulaire de 3 m x 6 m 
representee sur la Figure 7.4 
Reponse : P = 1235000./V 
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Huilc (d =r 0,S) 



Figure 7.3: 


* (Schaum 3.5 p.39) L’eau monte jusqu’au niveau E dans le conduit fixe au reservoir ABCD de la Figure 7.5. En 
negligeant le poids du reservoir et du conduit (a) determiner la force resultante agissant sur la surface AB qui a 2.50m de 
large, (b) determiner la force totale s’exercant sur le fond du reservoir et (c) comparer le poids total de l’eau avec le resultat 
(b) et expliquer la difference. 

Reponse : (a) P = 230 kN agissant ä 4.77m. sous le point O. (b) 839 kN ; (c) 298/clV. 

* (Schaum 3.7 p. 41) Le reservoir de la Figure 7.6contient de l’huile et de l’eau. Trouver 1a. force resultante agissant sur 
le cote ABC qui a 1.20m de large. 

* (Schaum 3.12 p. 44) Le cylindre de 2m de diametre de la Figure 7.7 pese 22.25 kN et a une longueur de 1.50m. 
Determiner les reactions en A et en B en ne tenant pas compte du frottement. 

Reponse : Reaction en A : 235401V dirigee vers la droite. Reaction en B : 37601V dirige vers le haut. 

* (Douglas 3.7 p.86) Une porte rectangulaire (Figure 7.8) d’epaisseur negligeable et de largeur b separe deux bassins 
ernplis du meine fluide densite p. On veut que la porte s’ouvre quand le niveau du bassin de gauche descencl en dessous d’une 
hauteur H par rapport ä Taxe de rotation. Le niveau du bassin de droite reste constant ä une hauteur y au dessus de Taxe. 
Exprimer le pds de la porte en fonction de p, g , b Y;y et H. 
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Figure 7.8: 


* (Schaum 4.5 p.61) Quelle est la fraction de volume d’un morceau de metal solide de densite 7.25 qui flotte ä la surface 
d’un recipient de mercure de densite 13.57 ? ( Figure 7.9) 

* (Schaum 4.8 p.62) De combien s’enfonce dans l’eau douce un tronc de 2;40?n de diametre et de 4.50m de long, de 
densite 0.425 ? ( Figure 7.10) 

* (Schaum 4.1 p.63) Un navire, dont les cotes sont verticaux ä la hauteur de la ligne de flottaison pese 4000 tonnes et a 
un tirant d’eau de 6.70m en eau salee ( p = 1025 kg/m 3 ). La decliarge de 200 tonnes d’eau de lest ramene le tirant ä 6.40m. 
Quel serait le tirant d’eau dans l’eau douce ? ( Figure 7.11) 

* (a)En negligeant l’epaisseur des parois du reservoir de la Figure 7.12, si le reservoir flotte dans la position indiquee, 
quel est son poids ? (b) si le reservoir est mainenu de Sorte que sapartie superieure soit ä 3 m au dessous du niveau de l’eau, 
quelle est la force qui s’y exerce (Figure 7.13)? 

* (Schaum 4.9 p.63) Une spliere d’aluminium pese dans l’air 0.258?V\ Son poids apparent dans l’eau n’est que de 0. 1001V. 
Montrer qu’elle est creuse et calculer son volume (rnasse volumique de l’aluminium 2580fcg/m. 3 ).(Rep.: V=58 cm 3 ). 

Application du theoreme de Bernoulli, pertes de Charge 

* Deversoir : un deversoir est un dispositif permettant de mesurer le debit qui s’ecoule dans un canal ä surface libre : on 
insere dans le courant une plaque que l’ecoulement doit donc surmonter pour s’ecouler. Dans cette plaque on dispose une 
fente de forme variable, ici triangulaire (cf.figure 7.14) et on mesure la hauteur h qui est fonction du debit. Exprimer le debit 
volumique en fonction de h et de 9. 

Reponse : Qv = tg(6/2) (j^h 5 ^ 2 ) 

* Conduite horizontale avec Venturi : un Circuit hydraulique horizontal a les dimensions indiquees sur la Figure 7.15. II 


47 






Figure 7.12: 



comporte un Venturi, une prise de pression statique et une prise de pression double. Ces divers appareils sont relies ä des 
manometres ä mercure. Les tubes de liaison aux manometres sont remplis d’eau. Masse volumique du mercure 13600fcg/m 3 . 
On donne : pression absolue en O : Pq = 1.5 bar. Pression atmosplierique : P a = 1 bar, z = Im. 

1 - Sacliant que la denivellation hy = 40mm calculer le debit volumique dans la conduite 

2 Calculer hs 

3 Calculer h d . 



Figure 7.15: Mesure d’un debit ä l’aide d’un Venturi 


* Dans le Circuit lrydraulique de la Figure 7.16 on donne: 
zi = 5m, Z 2 = zs = Z 4 , = Z 5 = 0, D 2 = -D 4 = 75 mm, D$ = 150 mm h = 10 cmdemer eure , y = Im. 
Coefficients de perte de Charge : £2 = 0.5 , £3 = 10. 



Calculer : 

1- Le debit volumique 

2- La pression effective ä l’aval des points 2, 3, 4. 
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* (Schaum 8.23 p.157) De l’huile circule du reservoir A ä travers 150?n de tuyau neuf de fonte asphaltee de 150mm de 
diametre jusqu’au point B de cöte 30m, comrne le montre la Figure 7.17. Quelle devra etre la pression en A pour que le 
debit de l’huile soit de 13 l/s (densite = 0.84 et v = 2.1 10 _6 m 2 /s )? Utiliser e = 0.12mm. 



Cole 86,0 m 


Figure 7.17: 


Machines hydrauliques 

* (Morel-Laborde p.172) Un laboratoire d’hydraulique dispose d’un banc oü un groupe electro-pompe P alimente, par 
l'intermediaire de conduite en acier galvanise, une turbine Pelton couplee ä un alternateur (cf. Figure 7.18) . La reserve 
d’eau se trouve sous la dalle du laboratoire. Pour le debit Q v = 2.5 l/s les caracteristiques de l’installation sont les suivantes 


• Groupe electro-pompe P : 

- rendement de la pompe rjp = 79% 

- rendement du moteur : 7j rn = 95% 

• Groupe turbine-alternateur T 

- rendement de la turbine : rjr = 82% 

- rendement de 1' alternateur : r/ a = 91% 

• Canalisations : diametre nominal 40 mm en acier galvanise e = 0.3mm 
DC\ = C 4 G = lm, C 2 C 3 = lm, BC 2 = 51m 

C\ et C 4 coudes 90° Rq = 2 DN 
C 2 et C 3 coudes 90° Rq = 1DN 

• Vanne ä opercule totalement ouverte 

• Sur le manometre M a 1 on lit une pression effective de 3.2 bar. 

1- Donner la puissance electrique fournie au groupe electro-pompe 

2- Donner la puissance electrique fournie par le groupe turbine-alternateur 

3- Quel est le rendement r]i de l’installation ? 

4- Quelle pression lira-t-on sur le manometre M a 2 
Pour l’eau :p = 1000 kg/m 3 ; v — 10 “ 6 m 2 /s. 

Reponse : 1 : P e \ = 1100 W ; 2- P e 2 = 438 W ; 3- r/- t = 40% ; 4- Pmo , 2 = 2.3 bar. 

• (Morel p. 158-159) Point de fonctionnement d’une turbine Pelton 

On veut se servir de l’eau du barrage de la Figure 7.19 pour produire de l’electricite. Pour cela on installe au bout de 
la conduite gravitaire une turbine Pelton couplee ä un alternateur dont les caracteristiques sont precisees dans le tableau 
suivant: 

H (m) 480 415 360 278 240 184 130 100 
Qv (l/s) 274 260 240 200 177 140 100 72 
r?(%) - 73 79.5 85 84 79 71 - 

Determiner graphiquement le point de fonctionnement de la turbine 
Reponse : H n = 333 m ; Qv = 228 l/s ; ryr = 82%. 
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Figure 7.18: 


T 1212 
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* (Morel p. 160) Point de fonctionnement d’une pompe avec pertes de Charge ä l’aspiration (Figure 7.20). La car- 
acteristique de la pompe est donnee Figure 7.21). 

1- Donner les points de fonctionnement de la pompe 2- Calculer la pression ä l’entree de la pompe 



Figure 7.20: Pompe avec perte de Charge a l’aspiration 


Continuite 

* Montrer que le champ de vitesse bidimensionnel defini en coordonnees cartesiennes par U x = kx ; U y = ky satisfait 
l’equation de continuite. 


Meme question en coordonnees cylindriques pour le champ de vitesse 1 D defini par : Ug = Ar + 


B 


La solution de l’ecoulement d’un fluide parfait incompressible autour d’un cylindre de rayon a s’ecrit : 

a 2 ' 

U r = Uncosd (l--r 


Ug = UosinO ^1 + - 

oü Uq designe la vitesse uniforme du fluide ä l’infini. Verifier que cette solution satisfait l’equation de continuite. 


On appelle source une singularite qui genere un debit de fluide. En coordonnees cylindriques une source produit un ecoulement 
radial dont les coordonnees sont : 

2nr 

Ug = U z = 0 


U r = 


Montrer que cette solution satisfait l’equation de continuite. 


* Soit le champ de vitesse defini par : 

U x = x 2 + y + z 
U y = x — xy + z 
U z = x + y — xz 

Verifier que ce champ de vitesse satisfait l’equation de continuite. 

* Soit un ecoulement permanent dont le champ de vitesse est defini par : 

U x = x + y 

U v = -x-y 
U z = 0 

Verifier que cet ecoulement se fait sans acceleration, c.a.d. que = 0. 

* Le champ de vitesse le long de la ligne centrale d’un convergent de longueur L est donne par : 

U x = 2t(l- 0.5 
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oü U x est la vitesse en m/s , t le temps en secondes mesure depuis le debut de la mise en mouvement du fluide et x la 
distance depuis l’entree du convergent. Trouver l’acceleration convective (Lagrangienne) et l’acceleration locale (Eulerienne) 
pour : t = 3s \ x = ^ et L = 0.8 m. 

* Le champ de vitesse bidimensionnel suivant satisfait il la condition de continuite ?: u = U {x 3 + xy 2 ) et v = U{y 3 + yx 2 ) 

* Meme question pour le champ suivant : a) u = 10 xt et v = —10 yt (p = Cste) b) u = U (y/S) 1 ^' et v = 0 (p = Cste) 

Theoreme de quantite de mouvement 


* Un trongon de tuyau coude se retrecit d’un diametre D i = 500mm en entree ä un diametre D 2 = 250 mm en sortie et 
procure une deviation de 45° ä l’ecoulement. Des mesures de pression indiquent P\ = 40 kPa en entree et P 2 = 23k Pa en 
sortie. Calculer l’intensite et la direction de la force resultante sur la conduite lorsque le fluide est de l’huile (p = 850kg/m 3 ) 
s’ecoulant avec un debit de 0.45m 3 /s. La conduite est horizontale. 

Reponse : R x = 5.454. 10 3 1V ; R y = 3.277. 10 3 iV ; angle =31°. 

* Un jet d’eau de diametre d = 50 mm s’ecoule ä la vitesse v = 4.9 m/s d’un orifice situe dans la paroi verticale d’un 
reservoir dont le niveau est maintenu a 1.5m au dessus du trou (Figure 7.22). Calculer la reaction du jet sur le reservoir a) 
lorsque le reservoir est immobile b) lorsque le reservoir se deplace en direction opposee au jet avec une vitesse constante de 
u = 1.2m/s. 



/ Volume de 
■ contröle 


V=4.9 m/s 


Figure 7.22: Effort sur un chariot 


Reponse : 47.14/V dans les 2 cas. 

* Un moteur-fusee consomme 1 kg de carburant pour 40 kg d’air traversant le moteur. La consommation de carburant 
est de 1.1kg /s quand l’avion se deplace dans de l’air au repos ä la vitesse de 200m/s. La vitesse des gaz issus du moteur 
ä la pression atmospherique est de 700m/s (vitesse relative par rapport au moteur). Calculer a)la poussee de l’engin, b) la 
puissance c) le rendement. 

Reponse : a) 22.77kN ; b) 4554/clU ; c) 44.7%. 

* La Figure 7.23 montre la coupe meridienne d’un ejecteur traverse par de l’air ( p = 1.2 kg/m 3 ) debouchant dans 
l’atmosphere ä travers une fente annulaire ä la vitesse de 30 m/s. Trouver la force s’exergant sur la conduite quand la pres- 
sion en A est inferieure de 2065Pa ä la pression de sortie. 

Reponse : 720.77V. 

Navier-Stokes 

* En coordonnees cartesiennes, ecrire les equations de Navier-Stokes d’un fluide circulant entre deux plans paralleles hori- 
zontaux avec les hypotheses suivantes (cf.figure 7.24): 

• ecoulement permanent 

• ecoulement etabli 

• ecoulement ID 

• fluide Newtonien incompresible 

* Idem en coordonnees cylindriques pour l’ecoulement dans un tube (cf.figure 7.25): 

* Idem en ecoulement tournant (cf.figure 7.26): 
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Figure 7.25: 



Figure 7.26: 
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* On considere l’ecoulcment laminaire d’un fluide de viscosite dynamique g = 0.9 Pa.s et de densite p = 1260 kg/m 3 entre 
deux parois d’extension infinie, separees de 10?7im et inclinees de 45° par rapport ä l’horizontale. La plaque superieure se 
deplace se deplace ä 1.5 m/s dans la direction opposee ä l’ecoulement. La plaque inferieure est au repos. Des jauges de pression 
installee en 2 points distants de 1 m verticalement enregistrent des pressions de 250kN.m~ 2 et 80 kN.m~ 2 . Determiner le 
profil de vitesse et de contrainte entre les plaques ainsi que le maximum de vitesse et la contrainte sur la plaque superieure. 



Figure 7.27: Ecoulement laminaire entre plaques paralleles. 


Reponse : U x = 566.4 y - 71.64.10 3 y 2 ; r y = 509.76 - 128.95. 10 3 y ; U max = 3.36 m/s ; r P = 0.78kN.m~ 2 . 

* On pompe de la glycerine de viscosite 0.9 Pa.s et de densite 1260 kg/m 3 dans un tube horizontal de 6.5m et de diametre 
d = 0.01m. Le debit vaut Q = 1.8 l/mn. Calculer le nornbre de Reynolds de l’ecoulement et determiner s’il est laminaire ou 
turbulent. Calculer la perte de clrarge. Calculer aussi le debit maximum possible en regime laminaire. 

Reponse : AP = 715 kN/m 2 ; Q cr it = 112Z/mn. 


* Fluides non-newtoniens (Guiffant p. 69) : certains fluides, dits de Bingham (cf. cours) comme des peintures ou meme 
de la boue ont la loi de comportement suivante : 

dv 

T <t c : -T- = 0 
dr 


T > T C : T = T C 


du 

■ a— = 0 
dr 


En decjä d’une certaine contrainte rq ces fluides se comportent donc comme des solides d’oü leur nom de fluide ä seuil. 
Trouver le profil de vitesse d’un fluide de ce type s’ecoulant dans une conduite cylindrique de rayon R en regime etabli et 
permanent sous l’effet d’un gradient de pression que l’on notera G. Exprimer le debit volumique Q en fonction de G, R , a 


et tc- 

Reponse : de r = 0 ä r = : U z (r) = Cste ; pour r > alors U z (r) = ^(r 2 i? 2 ) — ^-(r — R) 
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Formulaire 


Gradient d’un scalaire en coordonnees cylindriques ( r,9,z ): 


grad f 


Divergence d’un vecteur en coordonnees cartesiennes: divV = 

Divergence d’un vecteur en coordonnees cylindriques: divV = + ggf 

Continuite en cartesien : 

dp dpv x dpvy dpv z _ 
dt + dx + dy + dz 

Continuite en cylindrique: 

dp I 1 d{prv r ) , 1 d(pv e ) , d(pv z ) _ n 
dt ' r dr _r r de _r dz — u 

Equations de Navier Stokes en coordonnees cartesiennes (x, y, z) 


dl 

dr 

d±_ 


dz 


dv x dv x dv x 

~dT + Vx ~d^ + Vy ^ 


dVy 

dt 


dVy 

dx 

diu 


d'Vy 
+ v y~K — 

dy 


dv z 


dv x 

dz 

dvy 

dz 

dv z 


Equations de Navier Stokes en coordonnees cylindriques (r, d, z 


dp 

= -- 5 - + P9x 
dx 

dp 

dp 

= --z- + P9z 
dz 


P 


d 2 i 


d 2 i 


d 2 i 


dx 2 

d 2 v v 

dx 2 

d 2 v z 


dy 2 

d2v y | 
dy 2 

d 2 v z 

+ ~ — + 


dx 2 dy 


dz 2 

d 2 v y 

dz 2 

d 2 v z 

~dz 2 ~ 


dv r 

~dt 


dv r Vg dv r 


dv r 


r dr r d6 r z dz 


dp (Id dv r v r 

— äl + P9r + ß ( — ) - (2 


dr 


1 d 2 v r 2 dvg d 2 v r 


r dr y dr 


r 2 dd 2 r 2 dd dz 2 


'dvg dvg Vg dvg V r Vg dvg 

p[ ~dt +Vr ~d^ + v^e + ~ +Vz ~d. 7 


1 dp ( 1 d dvg Vg 1 d 2 Vg 2 dv r d 2 Vg 

' P9e + ß[ z — ( r —)~Z 5 + 


■de 


' dv z dv z vg dv 

P I — + V, 


diu 


r dr y dr ' r 2 ' r 2 dd 2 r 2 d9 dz 2 

1 d , dv z . 1 d 2 v z d 2 v ; 


dp 

r dd dz ) dz ^ PPz + ^ \ r dr y dr ’ r 2 dd 2 ' dz 2 


dt dr 

Contraintes visqueuses dans un fluide newtonien en coordonnees cartesiennes: 

_ 0 dv x 

T~XX rv 

OX 

T - 2u— 

Tyy ~ P dy 

diu 

T „ = 2 ß — 


(r-Z±) + 
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T~zr — T rz — 


dv z dv r 
dr dz _ 
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